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ABSTRACT 

In  this  report  a  simple  model  of  the  bubble  formed  by  detonating  a  limpet  mine  in 
contact  with  a  plane  surface  is  developed.  The  key  features  are  the  assumption  of  a 
hemispherical  bubble  and  the  account  taken  of  the  outflow  of  detonation  products 
from  the  bubble  through  the  hole  caused  by  the  detonation.  The  outflow  of  tike  gas  is 
assumed  to  be  one  dimensional.  The  equations  describing  the  time  variation  of  the 
bubble  radius  are  solved  for  various  values  of  tiie  physical  parameters  describing  the 
explosion.  The  model  displays  oscillatory  behaviour  which  is  damped  due  to  the  loss 
of  energy  via  the  outflowing  gas.  Under  some  circumstances,  the  outflow  of  gas  is  so 
large  that  only  one  oscillation  of  the  bubble  occurs. 
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On  the  dynamics  of  the  bubble  created 
upon  detonation  of  a  limpet  mine 


EXECUTIVE  SUMMARY 


When  an  explosion  occurs  underwater  there  are  two  principal  physical  phenomena 
that  follow.  A  shock  wave  is  transmitted  into  the  water  and  an  oscillating  bubble  is 
formed.  When  the  explosion  occurs  some  distance  from  the  target  the  damage  effects 
associated  with  these  two  phenomena  are  reasonably  well  understood.  When  the  deto¬ 
nation  occurs  in  very  close  proximity  to  the  target  the  damage  causing  effects  are  not  so 
well  understood. 

A  particular  underwater  weapon  that  detonates  in  very  close  proximity  to  a  target 
is  a  limpet  mine.  This  is  a  small  explosive  charge  that  is  detonated  in  contact  with  the 
hull  of  a  naval  target.  In  this  case,  not  only  does  a  shock  wave  form  and  propagate 
into  the  surrounding  ocean  and  through  the  target  but  the  explosion  itself  is  expected 
to  immediately  rupture  the  hull.  Subsequently  a  bubble  may  form  and  collapse  against 
the  bull.  Recent  studies  have  shown  that  when  an  explosion  bubble  collapses  against  the 
huU  of  a  naval  target  a  high  speed  liquid  jet  may  form  and  be  directed  towards  the  hull. 
This  effect  bears  a  degree  of  similarity  with  armour  defeating  shaped  charge  warheads. 

In  order  to  enhance  understanding  of  damage  causing  effects  associated  with  Hmpet 
mines  and  investigate  damage  mitigation  techniques,  a  series  of  simulated  detonations 
were  carried  out  against  the  decommissioned  destroyer  escort  HMAS  Derwent.  This 
series  formed  part  of  the  Ship  Survivability  Enhancement  Program  (SSEP)  conducted 
jointly  by  DSTO  and  the  RAN  in  November  1994. 

Part  of  the  instrumentation  of  the  limpet  mine  trial  series  of  detonations  was  high 
speed  underwater  photography,  in  order  to  gain  some  appreciation  of  the  possible  contri¬ 
bution  to  vessel  damage  by  the  collapse  of  the  bubble  formed  by  the  detonation.  Optimum 
deployment  of  the  photographic  equipment  required  an  estimate  of  the  maximum  bubble 
radius  and  period  of  oscillation  in  order  that  a  correct  field  of  view  and  filming  period 
were  employed. 

To  provide  these  estimates  an  analytical  model  of  the  motion  of  the  bubble  formed 
upon  detonation  of  a  limpet  mine  was  developed  and  is  described  in  this  report.  The 
model  takes  account  of  the  possibility  that  a  hole  may  be  formed  in  the  hull  upon  deto¬ 
nation  and  that  explosion  gases  may  escape  into  the  vessel.  The  calculations  show  that 
in  this  case  the  bubble  oscillations  are  damped  compared  to  the  behaviour  predicted  if 
no  hole  is  formed.  These  calculations  demonstrate  that  upper  estimates  of  the  maximum 
bubble  radius  and  bubble  period  may  be  made  using  the  model  and  assuming  that  no 
hole  is  formed  upon  detonation. 

Model  predictions  were  shown  to  be  in  reasonable  agreement  with  high  speed  photo¬ 
graphic  measurements  of  the  bubble  parameters.  It  is  concluded  that  the  model  provides 
a  useful  capability  to  predict  and  analyse  the  bubble  effect  of  limpet  mines  on  the  hull  of 
naval  platforms. 
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1.  Introduction 

A  limpet  mine  is  a  small  explosive  device  that  is  deployed  against  naval  targets.  It 
is  deployed  in  contact  with  the  target  and  magnets  are  usually  used  to  hold  the  mine  in 
place.  Limpet  mines  are  typically  deployed  in  a  clandestine  manner  by  divers.  During 
the  Ship  Survivability  Enhancement  Program  (SSEP)  conducted  by  DSTO  and  the  RAN 
during  November  1994,  a  series  of  simulated  limpet  mine  detonations  against  the  huU  of 
the  decommissioned  destroyer  escort  HMAS  Derwent  were  conducted  in  order  to  enhance 
understanding  of  the  damage  mechanisms  associated  with  these  weapons  and  possible 
techniques  for  damage  mitigation. 

Associated  with  an  underwater  detonation  are  two  principal  phenomena;  shock  wave 
formation  and  bubble  formation.  Typically  a  great  deal  of  emphasis  has  been  placed 
upon  the  damage  causing  potential  of  the  shock  wave  and  the  bubble  oscillation  has 
been  considered  in  the  context  of  hull  whipping,  which  is  primarily  a  far  field  bubble 
effect.  In  recent  years  bubble  collapse  in  close  proximity  to  targets  has  become  much 
better  understood,  with  the  phenomenon  of  high  speed  liquid  jet  formation  upon  bubble 
collapse  being  investigated  as  a  damage  causing  effect  possibly  comparable  to  the  shock 
wave.  In  view  of  this  progress  the  SSEP  limpet  mine  trial  series  was  instrumented  in  such 
a  manner  that  the  bubble  motion  event  could  be  visualised  using  high  speed  underwater 
photography  (Thornton  et  al.  1996).  In  order  to  determine  the  appropriate  field  of 
view  and  filming  time  requirements  a  simple  model  of  limpet  mine  bubble  oscillation  was 
developed  and  is  described  here. 

The  model  assumes  that  the  bubble  formed  upon  detonation  of  a  limpet  mine  is 
hemispherical  and  attached  to  the  vessel  hull  which  is  assumed  to  be  plane.  It  is  assumed 
that  a  fixed  proportion  of  the  explosion  energy  manifests  itself  in  the  bubble  motion. 
The  detonation  of  a  limpet  mine  will  almost  certainly  cause  the  formation  of  a  hole  in 
the  vessel  hull.  The  evidence  from  both  internal  and  external  high  speed  photography 
of  limpet  mine  detonation  events  during  the  SSEP  is  that  the  hole  is  formed  within  an 
interval  of  the  order  of  10  fxs  after  detonation  whereas  the  lifetime  of  the  bubble  is  of  the 
order  of  100  ms.  Some  account  of  the  effect  of  this  has  been  taken  in  the  model.  It  is 
assumed  that  a  hole  is  formed  immediately  upon  detonation  and  the  effect  is  that  bubble 
energy  is  lost  as  explosion  gases  are  transported  through  the  hole. 

The  equations  derived  for  bubble  radius  and  internal  pressure  are  solved  numerically 
to  gain  some  understanding  of  the  time  variation  of  the  limpet  mine  bubble  radius.  The 
radius  of  the  explosion  hole  is  varied  and  the  results  indicate  that  the  loss  of  energy  via 
gas  transport  through  the  hole  will  dampen  the  bubble  oscillation  and  as  the  hole  size  is 
increased  the  damping  increases  to  the  extent  the  bubble  does  not  rebound. 

Finally,  data  is  presented  for  the  variation  in  time  of  the  radius  of  the  bubble  formed 
following  one  of  the  limpet  mine  detonations.  In  this  example  the  adjacent  hull  compart¬ 
ment  was  flooded  and  some  comment  is  made  as  to  the  applicability  of  the  simple  model 
in  this  case.  This  single  experimental  data  set  is  in  acceptable  agreement  with  model 
calculations. 


2.  The  mathematical  model 

The  geometry  is  as  shown  in  figure  1.  The  centre  of  the  hemispherical  bubble  is  the 
origin  and  this  point  is  also  the  centre  of  the  hole  in  the  plane  boundary  caused  by  the 
detonation.  The  radius  of  this  hole  is  a  and  the  radius  of  the  hemispherical  bubble  is  R, 
which  is  time  dependent.  The  water  is  described  as  an  ideal  fluid  and  the  flow  irrotational 
and  hence  may  be  described  by  a  velocity  potential  4>-  Since  gravity  is  neglected  the  flow 
may  be  considered  as  due  to  a  source  of  strength  AttR^R  located  at  the  origin.  Hence  the 
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potential  is 


(1) 


where  r  is  the  position  vector  of  a  point  in  the  flow  field.  Note  that  this  velocity  potential 
satisfies  the  appropriate  rigid  boundary  condition  on  the  plane  boundary. 


Figure  1:  A  schematic  representation  of  the  bubble  formed  upon  detonation  of  a  limpet 
mine.  The  detonation  causes  a  hole  of  radius  a  to  be  formed  in  the  boundary. 


The  Bernoulli  equation  in  the  water  takes  the  form 

+  (2) 

where  p  is  the  pressure  in  the  fluid,  p^o  is  the  hydrostatic  pressure  at  the  point  of  detona¬ 
tion  and  p  is  the  water  density.  If  we  denote  by  pi  the  pressure  inside  the  bubble  which 
is  assumed  to  be  uniform,  then  evaluating  (2)  at  the  bubble  boundary  and  exploiting  (1) 
and  the  continuity  of  pressure  across  this  interface  yields  the  equation 

RR  +  l-R?  +  — =  0.  (3) 

Z  p 

In  the  case  where  there  is  no  hole  created  by  the  detonation,  so  that  no  explosion 
products  escape,  it  is  usually  assumed  that  the  products  behave  as  an  ideal  gas  undergoing 
adiabatic  expansions.  In  this  circumstance  p5  is  given  by 


P6  =  Po{Vo/Vy  =  po(Ro/R)^, 


where  V  denotes  the  volume,  the  subscript  0  denotes  initial  quantities  and  7  denotes  the 
ratio  of  specific  heats.  Discussion  of  the  determination  of  7,  po  and  Vq  is  postponed  until 
later  in  the  report.  It  is  noted  at  this  point  that  use  of  (4)  in  (3)  gives  the  Rayleigh-Plesset 
equation. 

In  order  to  determine  an  equation  for  the  time  variation  of  the  pressure  within  the 
bubble  it  is  necessary  to  consider  the  thermodynamics  of  the  bubble  contents,  taking 
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account  of  the  passage  of  gas  through  the  explosion  hole.  It  is  supposed  that  the  speed  of 
the  gas  flowing  through  the  hole  is  v  and  that  it  is  uniform  over  the  hole,  v  is  considered 
to  be  positive  when  the  gas  is  escaping  from  the  bubble.  The  detonation  products  are 
considered  to  be  ideal  so  the  equation  of  state  is 

pV  =  nRT,  (5) 

where  n  is  the  number  of  moles  of  gas,  R  is  the  universal  gas  constant  and  T  is  the 
temperature.  Since  heat  exchange  with  the  surrounds  is  neglected,  the  internal  energy  of 
the  gas  will  change  due  to  work  done  as  the  volume  changes  and  due  to  changes  in  the 
number  of  moles  of  gas  owing  to  outflow  and  inflow  through  the  hole. 

Under  these  circumstances  of  no  heat  exchange,  the  change  in  internal  energy  is 

dU  =  —p  dV  +  e  dn,  (6) 


with  e  the  internal  energy  per  mole  of  gas.  The  internal  energy  per  mole  of  gas  within 
the  bubble  is 


e  = 


RT 

7-1’ 


(7) 


which  is  the  same  as  the  internal  energy  per  mole  of  gas  escaping  or  entering,  since  it  is 
at  the  same  temperature.  The  volume  of  gas  escaping  per  unit  time  is 


so  the  number  of  moles  of  gas  escaping  per  unit  time  is 

dn  ■Kvno? 

li^  V~' 

The  internal  energy  of  the  gas  within  the  bubble  is 

nRT  pV 
U  =  - -  = 

7  —  1  7  —  1 

where  use  has  been  made  of  (5).  Hence  (6)  becomes 

— — (F  dp  pdV)  =  —pdV  +  edn. 
7-1 


Using  (6),  (7)  and  the  chain  rule  in  (9)  yields 


dp  TTvpa'^ 

V - h  7P  = - : — 

dV  V 


(8) 

(9) 


(10) 


(11) 


Noting  that  V  =  2tvR^/3  and  replacing  p  in  (11)  by  pf,  yields  the  equation  for  the  time 
variation  of  the  pressure  within  the  bubble  as 


dpj,  R  Sa^v 

IT  + 


(12) 


Notice  that  if  u  =  0  this  equation  may  be  immediately  integrated  to  give  (4). 

The  model  is  completed  with  specification  of  an  equation  giving  the  speed  with  which 
the  gas  escapes  from  the  hole.  An  idealised  situation  is  considered  in  order  to  expedite 
analysis.  It  is  supposed  that  the  gas  outflow  is  one  dimensional,  as  illustrated  in  figure  1. 
This  may  be  thought  of  as  assuming  flow  out  a  pipe  attached  to  the  hole.  The  outflow 
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is  assumed  isentropic  and  the  possibility  of  shock  wave  formation  is  neglected.  The 
assumption  of  isentropic  flow  permits  writing  the  pressure  as  a  function  of  the  density  as 


pp-^  =  K,  (13) 

where  /c  is  a  constant,  and  equal  to  that  value  that  describes  the  detonation  products  at 
t  =  0.  Since  the  flow  is  into  uniform  gas  at  atmospheric  pressure,  pa,  the  flow  is  a  simple 
wave  with  a  constant  Reimann  invariant.  For  outflowing  gas  this  condition  is  expressed 
as 

2  2 

- -c  -u  =  - -Co,  (14) 

7—1  7—1  ^  ' 

where  c  is  the  sound  speed  and  is  given  by 

=  ^p/p^  (15) 

and  Co  is  the  sound  speed  in  the  quiescent  fluid  in  the  far  fleld  where  the  gas  pressure  is 
equal  to  pa-  In  (14)  u  is  the  fluid  speed  in  the  outgoing  simple  wave. 

The  conditions  in  the  outflowing  gas  need  to  be  matched  to  the  conditions  within 
the  bubble  at  the  hole.  In  bubble  dynamics,  the  pressure  within  the  bubble  is  assumed 
uniform  and  the  flow  speed  of  the  gas  within  the  bubble  is  neglected.  Hence  it  is  assumed 
that  the  only  effect  that  outflow  through  the  hole  has  is  to  reduce  the  pressure  within 
the  bubble  in  accord  with  (12).  Thus  the  matching  of  conditions  is  achieved  by  assuming 
continuity  of  pressure  at  the  hole.  Then  c  in  (14)  may  be  evaluated  in  terms  of  pj,  and 
the  outflow  velocity  is  then  given  by  (14)  as 

» =  (16) 


where  use  has  been  made  of  (13)  to  eliminate  the  density  from  (15). 

To  proceed  to  numerical  solution  of  equations  (3)  and  (12)  it  is  helpful  to  introduce 
appropriate  scalings.  Distances  are  scaled  with  respect  to  Rm,  which  is  of  the  order 
of  the  maximum  bubble  radius.  Pressure  is  scaled  with  respect  to  pooj  the  hydrostatic 
pressure  at  the  detonation  point.  The  time  scale  is  Rmip/PooV^^  and  the  velocity  scale 
is  (p  oo! Scahng  with  respect  to  these  quantities  yields  the  equations  describing  the 
bubble  dynamics  as 


and 


RR  +  +  1  _  pj,  =  0, 

(17) 

dpb  „  R  ZcP'v 

(18) 

(19) 

non-dimensional.  It  should  be  noted  that  the  constant  k 

appearing  in  (19)  is  here  non-dimensionalised  in  accord  with  the  scaling  given  above. 


3.  Initial  conditions 

The  initial  conditions  utilised  for  studies  of  the  bubble  dynamics  are  that  the  bubble 
initially  has  a  non-dimensional  radius  Rq  and  that,  in  the  absence  of  gas  loss  through 
an  explosion  hole,  would  expand  to  a  maximum  non-dimensional  radius  of  one.  Given 
the  initial  bubble  radius  the  initial  bubble  pressure  may  also  be  determined  in  a  manner 
described  in  what  follows.  The  initial  radial  velocity  of  the  bubble  surface  is  taken  to  be 
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zero.  With  such  assumptions,  at  both  maximum  and  minimum  radius  the  kinetic  energy 
of  the  water  is  equal  to  zero.  Hence  the  total  energy  is  equal  to  the  sum  of  potential  and 
internal  energy.  The  potential  energy  is  the  work  done  against  the  hydrostatic  pressure 
in  forming  a  bubble  of  radius  R  and  is  equal  to 

2 

Poo) 

and  the  internal  energy  is  the  work  done  in  compressing  the  bubble  contents  adiabatically 
from  infinite  volume  to  the  bubble  V  and  is 


The  statement  of  energy  conservation  at  maximum  and  minimum  volume  is  then 

+  =  (20) 

3  3  7  —  1 

where  Eq  is  the  energy  in  the  bubble  motion.  In  this  expression,  the  quantity  pqRq^ 
may  be  eliminated  in  favour  of  k  given  in  (13).  In  order  to  perform  calculations  for 
practical  problems,  this  is  the  strategy  followed,  with  k  determined  from  experimental 
observations. 

To  fix  ideas,  consider  that  the  explosive  is  TNT.  Experimental  data  recorded  in 
Taylor  (1950)  indicates  that  the  ratio  of  specific  heats  may  be  taken  as 

7  =  5/4  (21) 


and  that 

pp-'T  =  K  =  1.386  X  10^  (22) 

where  p  and  p  are  respectively  measured  in  kgm“^  s~^  and  kgm~^.  The  energy  yield  per 
kilogram  of  TNT  is  denoted  by  £  and  is  equal  to  3.683  x  10®  Jkg“^  so  supposing  that 
a  proportion  p  of  the  total  energy  yield  of  the  explosion  manifests  itself  in  the  bubble 
motion,  the  energy  balance  at  maximum  and  minimum  radius  is 

\^R%oo  +  (^)i-'>-^ie-®(^-i)  =  pw£,  (23) 

3  3  7  —  1 

where  w  is  the  charge  mass  in  kilograms.  In  the  calculations  presented  later  p  is  taken 
to  be  equal  to  1/2. 

It  is  convenient  at  this  stage  to  introduce  a  distance  scale 


/3ptyi5\i/3 

V  27rpoo  ^ 


(24) 


and  if  R  =  ^R*  then  (23)  becomes 


with 


R*^  +  =  1,  (25) 

'=(7^)(f‘7'^J’~''  <26) 


which  is  a  small  number.  The  maximum  and  minimum  radii  of  the  bubble  are  then  given 
by 


Rmax  — 

■Rmin  ~  -^2^’ 


(27) 
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where  and  R2  are  the  roots  of  (25)  which  are  given  with  good  accuracy  as 

=  1  -  e/3  -  7e2/36, 

R^  -  e^/2 

Having  determined  the  maximum  and  minimum  radii,  it  remains  to  determine  the  non- 
dimensional  value  ofpo>  the  initial  pressure  within  the  bubble.  Using  (22)  and  (28),  the 
dimensional  value  of  po  is 

PO  =  K,p'l^{p£)~^  R2~^  ,  (29) 

and  the  non-dimensional  value  of  this,  also  known  as  the  strength  parameter  a  (Best  &; 
Kucera  1992),  is 

(30) 

The  initial  non-dimensional  bubble  radius  is 

Rq  —  RjDin/ Rmdx.  — R2I  ^1-  (31) 


4.  Example  calculations 

In  order  to  integrate  equations  (17)-(19)  a  fourth  order  Runge-Kutta  method  has 
been  employed  using  a  variable  time  step.  The  time  step  is  given  by 

A 

At  =  - - - , 

R2/2  +  P&+1 

where  is  a  constant  here  taken  to  be  equal  to  0.01.  This  choice  of  time  step  has 
been  introduced  by  Best  &  Kucera  (1992)  in  order  to  perform  boundary  integral  method 
calculations  of  non-spherical  bubble  motion.  This  choice  of  time  step  ensures  that  small 
steps  are  taken  when  the  speed  of  the  bubble  boundary  is  large  or  when  the  bubble 
contents  are  highly  compressed  near  to  rebound.  At  other  times  larger  steps  are  taken 
to  minimise  the  computational  resources  utilised. 

To  illustrate  the  features  of  the  model  a  number  of  calculations  are  presented  for 
various  charge  masses  and  various  explosion  hole  diameters.  This  model  provides  no 
estimate  of  the  size  of  the  explosion  hole  expected  to  be  created  by  the  detonation  and 
this  quantity  is  simply  input  to  the  model.  The  results  presented  here  are  thus  only 
indicative  of  the  behaviour  that  might  be  expected  in  a  real  detonation  scenario.  Figures 
2-4  illustrate  the  variation  of  bubble  radius  with  time  for  charges  of  mass  10  kg,  5  kg  and 
1  kg,  these  being  values  that  were  considered  for  the  SSEP  limpet  mine  trial.  The  depth 
of  detonation  is  3.5  m  which  corresponds  to  deployment  of  the  mine  near  the  keel  of  a 
typical  target  warship. 

Given  the  strong  assumptions  made  in  derivation  of  the  model  equations,  particularly 
concerning  the  flow  of  the  explosion  gases  through  the  hole,  these  results  are  best  used 
to  determine  qualitatively  the  trends  in  bubble  behaviour  with  varying  charge  mass  and 
hole  radius.  The  trend  with  respect  to  increasing  hole  size  is  illustrated  by  all  examples 
and  figure  2  is  as  good  as  any.  For  the  case  where  no  explosion  hole  is  formed  there 
is  no  mechanism  in  the  model  whereby  energy  is  lost  from  the  bubble  motion.  Thus  an 
oscillatory  motion  is  predicted  as  illustrated  in  this  figure.  In  reality  any  non-uniformity  in 
the  flow  field  surrounding  the  bubble  will  cause  the  bubble  to  deform  from  hemispherical 
shape  and  a  liquid  jet  will  most  likely  form  upon  the  collapse.  Such  sources  of  flow 
asymmetry  include  buoyancy  and  the  flow  field  geometry  itself.  It  should  be  noted  that 
the  expansion  phase  of  a  bubble  motion  is  stable  and  any  deviations  from  hemispherical 
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shape  are  expected  to  occur  during  the  later  stages  of  the  bubble  collapse.  Thus  in  the 
event  that  a  hole  is  not  formed  by  the  detonation  this  model  is  expected  to  well  predict 
the  motion  of  the  bubble  for  most  of  its  first  oscillation  and  particularly  its  maximum 
radius  and  period. 

As  the  hole  caused  by  the  detonation  increases  in  radius  the  oscillation  of  the  bubble 
becomes  damped.  In  figure  2  the  hole  of  radius  0.05  m  causes  the  first  maximum  to  be 
less  than  in  the  case  of  no  hole.  The  bubble  rebounds  in  this  case  with  the  rate  of  collapse 
near  the  minimum  so  great  that  despite  the  outflow  of  gas  through  the  hole  the  pressure 
builds  to  an  extent  that  it  cannot  be  relieved  by  gas  outflow  and  the  bubble  rebounds. 
The  second  maximum  is  damped  with  respect  to  the  first  as  is  the  third  maximum  with 
respect  to  the  second.  As  the  bubble  collapses  a  third  time  the  remaining  gas  within  the 
bubble  flows  out  the  hole  and  the  bubble  collapses  completely. 


T  ime  ( s ) 

Figure  2:  The  bubble  radius  as  a  function  of  time  following  detonation  of  10  kg  TNT  at 
a  depth  of  3.5  m.  The  curves  correspond  to  an  explosion  hole  of  radius:  (a)  0.0  m,  (b) 
0.05  m,  (c)  0.1  m,  (d)  0.25  m. 


The  final  two  examples  illustrate  the  case  where  the  hole  is  so  large  that  all  bubble 
gases  flow  out  during  the  first  collapse  and  the  bubble  does  not  rebound.  It  is  interesting 
to  note  the  shape  of  the  radius-time  curve  for  the  case  a  —  0.25  m.  All  other  oscillation 
periods  of  the  bubble  are  qualitatively  symmetrical  about  the  maximum  but  in  this  case 
the  curve  is  noticeably  skewed. 

The  examples  for  smaller  charge  masses  illustrated  in  figures  3  and  4  show  similar 
features.  As  the  hole  size  increases  the  bubble  motion  is  damped  to  the  extent  that  the 
bubble  eventually  collapses  completely.  Notice  that  in  the  case  where  tu  =  1  kg  and 
a  =  0.025  m  the  third  and  final  bubble  period  demonstrates  significant  asymmetry  in  the 
radius-time  curve. 
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Figure  4-  The  bubble  radius  as  a  function  of  time  following  detonation  of  1  kg  TNT  at 
a  depth  of  3.5  m.  The  curves  correspond  to  an  explosion  hole  of  radius:  (a)  0.0  m,  (b) 
0.025  m,  (c)  0.05  m,  (d)  0.1  m. 
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This  model  has  been  formulated  making  some  strong  assumptions  as  to  the  flow  of 
gas  through  the  explosion  bubble  hole.  Indeed,  if  the  variation  of  escaping  gas  velocity 
with  time  is  plotted  as  in  figure  5,  for  the  case  where  u;  =  10  kg  and  a  =  0.05  m,  then 
the  high  values  near  bubble  minimum  suggest  that  shock  waves  may  be  formed.  Note 
also  that  when  the  bubble  is  near  its  maximum  radius  the  pressure  within  the  bubble  is 
so  low  that  the  model  predicts  gas  flow  into  it.  The  possibility  of  shock  formation  in  the 
outflowing  gas  could  be  incorporated  into  the  model,  but  given  the  lack  of  information 
regarding  key  aspects  of  the  limpet  mine  detonation  event,  such  as  the  formation  of  the 
hole  and  the  proportion  of  the  explosion  energy  yield  that  manifests  itself  in  the  bubble 
motion,  such  effort  seems  unwarranted.  To  proceed  detailed  hydrocode  modelling  would 
appear  to  be  the  most  appropriate  course. 


T  ime  (  s  ) 


Figure  5:  The  escaping  gas  velocity  as  a  function  of  time  for  the  example  of  detonation 
of  10  kg  TNT  at  a  depth  of  3.5  m.  The  explosion  hole  is  of  radius  0.05  m. 


Table  1:  The  predicted  upper  bounds  on  the  maximum  bubble  radius  and  first  oscillation 
period  for  a  range  of  charge  masses.  All  detonations  are  assumed  to  take  place  at  a  depth 
of  3.5  m. 


Charge  mass  (kg) 

Maximum  radius  (m) 

Oscillation  period  (s) 

10 

3.79 

0.63 

5 

3.01 

0.50 

1 

1.76 

0.29 

The  primary  aim  of  the  calculations  undertaken  here  was  to  make  estimates  of  the 
period  and  maximum  radius  of  the  bubble  formed  upon  detonation  of  a  limpet  mine. 
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In.  this  context  the  model  case  where  no  hole  is  formed  by  the  detonation  is  expected  to 
provide  upper  bounds  on  these  quantities,  with  the  calculations  performed  here  indicating 
that  the  effect  of  hole  formation  is  to  dampen  the  bubble  oscillation.  Hence  for  the 
charge  masses  considered  here  and  detonated  at  a  depth  of  3.5  m  table  1  shows  the  upper 
estimates  of  the  period  and  maximum  radius. 

During  the  SSEP  limpet  mine  trial  high  speed  underwater  photography  was  per¬ 
formed  for  a  number  of  detonations  (Thornton  et  al.  1996).  The  photographic  record 
of  a  single  shot  was  amenable  to  detailed  analysis  and  yielded  the  radius-time  history 
shown  in  figure  6.  This  example  is  for  2  kg  of  PE4  detonated  at  a  depth  of  2  m.  It  is 
emphasised  that  in  this  case  the  internal  ship  compartment  adjacent  to  the  detonation 
point  was  completely  flooded  with  sea  water.  Measurements  made  after  the  event  showed 
a  hole  of  approximate  diameter  275  mm  compared  with  the  charge  diameter  of  170  mm. 
Further,  the  hole  showed  a  light  secondary  indentation,  indicating  that  the  hole  size  was 
not  significantly  increased  by  a  flow  of  water  or  gas  through  it  during  the  bubble  motion. 
Hence  it  may  be  assumed  that  the  hole  was  formed  essentially  upon  detonation. 

This  scenario  is  not  considered  to  be  appropriately  described  by  the  model  developed 
here.  However,  if  the  case  is  considered  where  a  limpet  mine  is  detonated  adjacent  to 
a  plane  boundary  surrounded  on  both  sides  by  water  then  it  may  be  argued  that  the 
bubble  oscillation  parameters  may  be  predicted  assuming  motion  of  a  bubble  in  the  free 
field.  If  it  is  considered  that  the  shock  formed  upon  detonation  propagates  only  into 
the  half  space  on  which  the  charge  is  located  and  that  the  energy  that  would  otherwise 
manifest  itself  in  a  shock  propagating  into  the  other  half  space  is  essentially  expended  in 
creating  the  explosion  hole,  then  it  may  be  deduced  that  the  remaining  energy  manifests 
itself  in  the  bubble.  In  accordance  with  the  evidence  supporting  the  proposition  that  the 
hole  is  formed  immediately  upon  detonation,  there  is  no  preferential  direction  in  which 
the  explosion  gases  will  expand.  Thus  the  gases  will  expand  equally  on  either  side  of  the 
boundary  and  the  bubble  so  formed  will  be  equivalent  to  that  formed  if  the  detonation 
occurred  in  the  free  field.  In  this  case  the  statement  of  energy  conservation  is  slightly 
different  owing  to  the  spherical  shape  of  the  bubble  with  the  potential  energy  given  by 

^ttR^Poo, 


and  the  internal  energy  is 


87—1 


The  analysis  foUows  through  as 
length  scale  j3  is  given  by 


for  the  hemispherical  bubble  with  the  exception  that  the 


(33) 


The  non-dimensional  maximum  and  minimum  radii  are  determined  as  the  roots  of  (25) 
and  the  expression  for  the  non-dimensional  parameter  e  is  exactly  (26).  The  initial 
bubble  pressure  and  strength  parameter  are  given  by  (29)  and  (30)  and  the  initial  non- 
dimensional  bubble  radius  by  (31). 

In  the  absence  of  better  data  it  is  assumed  that  the  performance  of  PE4  underwater 
may  be  described  by  the  same  parameters  as  TNT.  Evaluating  the  maximum  radius  and 
bubble  period  using  this  model  yields  iJmax  =  1-83  m  and  T  =  0.32  s.  These  values  agree 
favourably  with  those  observed  as  evident  from  inspection  of  figure  6.  Note,  however, 
that  the  correspondence  is  particularly  good  during  the  early  phase  of  the  expansion  but 
that  the  experimental  data  fall  away  from  the  theoretical  curve  as  the  bubble  collapses. 
It  was  observed  that  the  bubble  in  this  case  did  not  rebound  so  this  trend  in  the  bubble 
radius  is  likely  to  be  due  to  the  outflow  of  explosion  bubble  gases  into  the  target.  Note 
that  the  asymmetry  of  the  radius  time  curve  bears  qualitative  similarity  to  some  of 
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the  calculated  radius  time  curves  for  hemispherical  bubbles  for  those  oscillation  periods 
where  the  bubble  collapses  completely.  This  single  example  provides  some  evidence  that 
significant  limpet  mine  bubble  parameters  may  be  approximately  determined  using  the 
elementary  analysis  presented  here. 


a 


;:3 

K 


T  ime  (  s  ) 


Figure  6:  The  bubble  radius  as  a  function  of  time  as  determined  from  a  high  speed  pho¬ 
tographic  record  of  the  detonation  of  2  kg  of  PE4  at  a  depth  of  2  m.  In  this  example  the 
internal  ship  compartment  adjacent  to  the  detonation  point  was  flooded.  The  experimen¬ 
tal  data  are  denoted  by  a  cross  and  the  solid  curve  is  calculated  assuming  motion  of  a 
spherical  bubble. 


5.  Conclusions 

In  this  report  an  analytical  model  has  been  developed  that  describes  the  motion  of  the 
bubble  formed  upon  detonation  of  a  limpet  mine  adjacent  to  a  plane  boundary.  The  model 
allows  predictions  to  be  made  of  the  maximum  bubble  radius  and  period.  A  significant 
factor  that  affects  the  behaviour  of  the  bubble  is  the  radius  of  the  hole  presumed  formed 
upon  detonation.  The  hole  allows  explosion  bubble  contents  to  vent  into  the  vessel 
hull  and  thus  dampens  the  bubble  oscillation.  The  results  of  calculations  indicate  that 
for  a  small  hole  a  number  of  bubble  oscillations  may  occur,  with  the  maximum  radius 
decreasing  for  each  subsequent  oscillation  period.  As  the  hole  radius  increases  a  stage 
is  reached  where  ail  the  bubble  contents  are  vented  on  the  first  collapse  and  the  bubble 
does  not  rebound. 

A  second  significant  factor  aflFecting  bubble  behaviour  is  the  manner  in  which  the 
energy  yield  of  the  explosion  manifests  itself  in  propagating  a  shock  wave  into  the  water, 
cutting  a  hole  in  the  vessel  hull  and  in  bubble  motion.  It  has  been  assumed  here  that  half 
the  available  energy  manifests  itself  in  bubble  motion.  In  the  absence  of  any  data  this 
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may  be  a  reasonable  first  approximation  but  is  an  issue  requiring  further  investigation  if 
more  detailed  calculations  are  required. 

Finally,  the  primary  aim  of  this  modelling  exercise  was  to  provide  upper  bounds  on 
the  maximum  bubble  radius  and  period  for  detonation  of  a  limpet  mine  against  the  hull 
of  a  naval  target.  It  is  expected  that  use  of  equations  (17)-(19)  with  the  hole  radius  set 
equal  to  zero  will  provide  appropriate  upper  bounds.  The  effect  of  a  non-zero  hole  radius 
will  be  to  dampen  the  motion  compared  to  these  values. 
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